Diffusion of a sphere in a dilute solution of polymer coils 
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We calculate the short time and the long time diffusion coefficient of a spherical tracer particle in a polymer 
solution in the low density limit by solving the Smoluchowski equation for a two-particle system and applying 
a generalized Einstein relation (fluctuation dissipation theorem). The tracer particle as well as the polymer coils 
are idealized as hard spheres with a no-slip boundary condition for the solvent but the hydrodynamic radius of 
the polymer coils is allowed to be smaller than the direct-interaction radius. We take hydrodynamic interactions 
up to 1 1th order in the particle distance into account. For the limit of small polymers, the expected generalized 
Stokes-Einstein relation is found. The long time diffusion coefficient also roughly obeys the generalized Stokes- 
Einstein relation for larger polymers whereas the short time coefficient does not. We find good qualitative and 
quantitative agreement to experiments. 



I. INTRODUCTION 

Transport properties of Brownian particles in suspensions 
are of great interest for all technological applications involv- 
ing complex fluids such as food technology or oil recovery 
and they have been studied extensively experimentally (see, 
e.g., SH) and theoretically (see, e.g., SHUSH). The 
diffusion constant D s of a spherical tracer particle with radius 
R s in a simple solvent is to a good approximation given by 
the Stokes-Einstein relation 



pendent diffusion coefficient D s (t) via 12111 

D s {t)t= i((r s (t)-r s (i = 0)) 2 ), (2) 
6 

with the particle position r s (t) at time t. (•} indicates the equi- 
librium ensemble average. The short and the long time limit of 
D 8 (t) are called the short and long time diffusion coefficients 



D s „ = limDJt) and 



Dt 



D s = 



k B T 
6 7r R s i] Q ' 



lim D s (t), 

t—>oo 



(3) 
(4) 



(1) 



with the solvent viscosity 770 and the thermal energy fcg T. As 
demonstrated, e.g., for the case of a tracer sphere in solution of 
polymers S SHU EL me naive approach to replace the pure 
solvent viscosity with the macroscopic shear viscosity T] maao 
of the polymer solution (as measured in a viscosimeter) in 
general fails. The polymer solution in the vicinity of the mov- 
ing sphere is not homogeneous. Even in equilibrium one ob- 
serves depletion layers or density oscillations (depending on 
the interaction potentials between the polymers and between 
the polymer and the particle). In the vicinity of a moving par- 
ticle, the flowing solvent rearranges the polymers leading to 
an enhanced polymer density in-front and a reduced polymer 
density behind the particle fl2ll . which leads to an enhanced 
friction ifl^fl^fTstl . and to long-ranged solvent mediated ef- 
fective interactions HH [l7tl . The time scale for the build-up 
of these inhomogeneities in the solution in the vicinity of the 
moving particle is given by the diffusivity of the polymers and 
the particle size. For most systems this time scale is well sep- 
arated from the corresponding microscopic time scale of the 
solvent lfl8l [l9i I20I1 . but rather close to the time scale of the 
particle diffusion (given by the particle size and its diffusion 
constant in the pure solvent). 

For such systems, the mean square displacement is not lin- 
ear in time and for the tracer particle one defines a time de- 



respectively. The diffusion coefficient D s (t) is related to the 
time-dependent mobility coefficient fj, s (t) via the fluctuation- 
dissipation theorem, i.e., the generalized Einstein relation, 



(5) 
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with j3 = 1 / ( ks T) . Hs(t) is the linear response mobility de- 
fined by the ratio of the average velocity of the particle and 
a small and constant external force F ext that starts to act on 
the particle at t = 0. Therefore, at time t = 0, the distri- 
bution of polymers around the sphere is still in equilibrium, 
and the short time mobility = limt_>o Ms (i) = P D s s is 
solely determined by the hydrodynamic forces on the tracer. 
Once the sphere is in motion, the distribution of polymers 
in the vicinity of the tracer becomes anisotropic: it is more 
likely to find a polymer in front of the sphere than behind it, 
which reduces the mobility of the sphere. After a sufficiently 
long time, the polymer distribution becomes stationary and 
the velocity is related to the force via the long time mobility 
/4 = lim^oc fj, s (t) = (3 D l s . In general (i l s < /if. 

Most of the theoretical work focuses on the semi-dilute 
regime with effective or mean field models [22, 23. 24. 25]. 
see rf26ll for a summary. The depletion of polymers near the 
surface of the colloid was considered in Ref. [27]. Recently, 
a model based on the reduced viscosity of the polymer solu- 
tion near the colloid due to depletion was introduced l28ll29ll . 
However, up to now, the distinction between short and long 
time diffusion coefficient has been considered only for sus- 
pensions of hard spheres S US HH HI HI . One main dis- 
tinction between hard spheres and polymers is, that for hard 
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FIG. 1: A tracer sphere with radius R s is suspended in a dilute so- 
lution of polymer coils with radii of gyration R g . In a model with 
hard-sphere interactions, the centers of mass of the polymer coils 
cannot pass the dashed surface with radius R = R g + R s . 



direct 

interactions: 

hydrodynamic 
interactions: 





FIG. 2: We model the polymer coils as hard spheres of radius R g 
concerning the interactions with the tracer sphere sphere and as a 
solid sphere of radius Rh and a no-slip boundary condition on the 
surface concerning the interactions with the solvent. 



spheres the hydrodynamic radius is equal to the particles ra- 
dius, while for polymers, the effective hydrodynamic radius 
is in general smaller than the radius of gyration. A different 
hydrodynamic and sphere radius has been used as a model 
for charged colloids l34ll . but only the case of equal tracer 
and bath particles was considered. In this paper, we focus 
on the regime of a dilute polymer solution but we distinguish 
between the short and long time diffusivity. We idealize the 
tracer particle and the polymer coils as hard spheres concern- 
ing both, the direct and hydrodynamic interactions, but with 
a hydrodynamic radius which can be smaller than the interac- 
tion radius in the case of the polymers. 

In the following section, we present our model system and 
in Sec. [II] we calculate the short and long time diffusion co- 
efficients from the corresponding Smoluchowski (or Fokker- 
Planck) equation. The results are compared to experimental 
values in Sec. Mand we conclude in Sec. [VI] 



II. MODEL 

We model the tracer particle as well as the polymer coils 
as spherical overdamped Brownian particles with radii R s 
and R g , respectively, as shown in Fig. Q] The bare diffusion 
coefficients of the sphere and the polymers in the pure sol- 
vent are D a s and D p (Indices s and p will denote the sphere 
and the polymer, respectively, throughout the paper.), respec- 



tively, which are calculated via the Stokes-Einstein relation in 
the pure solvent with viscosity r/ . The hydrodynamic radius 
which enters the Stokes-Einstein relation is equal to R s for the 
tracer and Rh < R g for the polymers, see Fig. [2] We idealize 
the direct interaction between the tracer and the polymers as a 
hard-sphere interaction 



V(r) = 



for r> R a + R g 
oo for r < R s + R„ 



(6) 



r = |r| = |r s — r p | denotes the distance of the centers of the 
tracer and a polymer. As a consequence, the center of mass of 
the polymers can approach the center of the tracer only up to 
a distance R = R s + R g . In the dilute limit, we neglect the 
mutual interactions of the polymer particles, which reduces 
the problem effectively to a two-particles system of one tracer 
sphere and one polymer coil. 

The dynamics of the system is described by the Smolu- 
chowski equation for the probability density P(r s ,r p ) for 
finding the tracer sphere at position r s and the polymer coil 
at r p . The Smoluchowski equation is a continuity equation 
and with the corresponding probability currents j s i p we can 



define the velocity operators such that j 



d 



v, 

V, 



s/p - 



P: 
v„P 



D 



(3 [V P V] + V f 



P 



, (7) 



with the external force F ext . V s and V p denote the gradient 
with respect to the position of the sphere and the polymer, re- 
spectively. The components of the symmetric diffusivity ma- 
trix 



D 



ss 

sp 



(8) 



have the form 

D 
D 



i>i> 



D 



ps 



D° s [D r ss (r)P+D e ss (r)(l-P)] (9a) 
D p [D r pp (r)P + D e pp (r)(l-P)} (9b) 
D sp = D° s [D; s (r)P + D e ps (r)(l - P)] . (9c) 



The projector P is given by rf, with r = r/r. The r- 
dependent coefficients D r as , D s ss , D pp , D pp , D p8 , and D ps can 
be expanded in a power series in r -1 . We use the coefficients 
for spheres with no-slip hydrodynamic boundary conditions 
on their surfaces up to order r~ n in the distance according to 
Ref . ll35ll . We quote the coefficients for two unequal spheres 
of radii R s and Rh for the tracer and the polymer, respectively 
in appendixlAl We therefore assume that the polymer interacts 
with the solvent like a solid sphere with radius Rh, see Fig. 
[2] For the case of hard sphere suspensions, also lubrication 
forces have been taken into account 15[]. In the case consid- 
ered here, the radius of gyration R g is always larger than Rh, 
i.e., the hydrodynamic ally interacting spheres never come into 
contact and the far field expansion converges well. The diffu- 
sion matrices (O are valid on the Brownian time scale and for 
small Reynolds numbers l2~lll . 
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Eq. (0 is translationally invariant since D depends only on 
r = r s — r p . As a consequence, P is a function of r only, and 
the hard interaction potential V in Eq. (0 can be translated 
into a no-flux boundary condition on a sphere of radius R 



(is 



at 



R. 



(10) 



In thermal equilibrium (for which F ext = is a necessary 
condition) detailed balance holds and all components of the 
probability currents j s / p are zero. The equilibrium distribu- 
tion and therefore also the initial condition for the dynamical 
problem Eq. (0 is therefore given by 

P(r, t) = P eq {v) = p 9(|r| - R) at t = 0, (11) 

with the average number density of polymer molecules p. Far 
from the tracer sphere, the polymer distribution should be un- 
affected by the presence of the sphere and therefore equal to 
the corresponding equilibrium distribution, which yields the 
boundary condition 



P(r,t) 



P 



for 



(12) 



In the long time limit P(r, t) reaches a steady state 
P(r,t) — > P°°(r) for t — > oo, which is given by the solu- 
tion of the stationary version of Eq. (0 (for r > R) 

= V • (D ss + D pp ) • VP 00 - V • D ss • P°° F ext 

- V • [D ps • {2VP°° -/3P°°F ext }] . (17) 

Since we are interested in the linear response regime, we ex- 
pand P°° (r) in powers of F ext up to linear order, i.e., we seek 
a solution of the form lull 



(18) 



P°°(r) = P e «(r) [1 + pA(r)i ■ F ext ] 



Inserting Eq. (IT8b into Eq. ( TPTl i and keeping only terms linear 
in F ext yields a second order linear differential equation for 
the coefficient A(r). Since we expand the mobility matrix 
D(r) in Eq. (0 in a power series in 1/r up to order 0(r~ n ), 
we choose the following ansatz for A(r) 



(19) 



In the linear response regime, the average velocity of the 
tracer particle is given by ((-) F (t) denotes the time dependent 
non-equilibrium average) 

(v s ) F (t)= /j.(r,t) d 3 r = Jv s P(r,t) d 3 r = » a (t) F ext , 

(13) 

such that we can calculate the short and long term diffusion 
coefficients from the solution of the Smoluchowski equation 
©: once P(r, t) is known, Eq. (TUt yields the mobility coeffi- 
cient /i s (t), from which we calculate the diffusion coefficients 
through the generalized Einstein relation Eq. (0. Since D s (t) 
has to be finite for t — ► 0, the integration constant which ap- 
pears when solving Eq. Q for D s has to be zero. 

Inserting the expression for v s P from Eq. (0 into Eq. ( fT3l 
we can decompose the velocity into three components 

(v s ) F (t) = (3 (D ss > ■ F** + (vi) F (t) + (vf ) F (t), (14) 



with 



-/3(D SS -VF-D sp - VV) F (t), (15) 



(vf) F (t) = 

vf r ) F W = -/3(D ss -VlnP-D sp - VlnP) F (t). 



(16) 

Since P only depends on r = r s — r p , we have replaced the 
gradients with respect to the positions of the tracer and the 
polymer with V = V r , = — V rp . and v^ r are the result 
of the direct interactions and the Brownian force, respectively 

El. 

In the short time limit, i.e., at t = 0, P(r, t) is given by the 
initial condition, i.e., by P eq (r). By symmetry, both (vf) and 
(vf r ) are zero in this limit. 



which turns this differential equation into an algebraic equa- 
tion for the coefficients c;. It can be solved for the c^2 in 
terms of c 2 . Note that C\ = 0. c 2 finally is determined by 
the boundary condition ( TTOb at short distances, which reads in 
terms of A(r), 



dA(r) 



dr 



D r ss (R) - DUR) 



-R D r ss (R) 



^D; p (R) - 2Dl p (R) 



(20) 



We solve the above equations using the computer algebra sys- 
tem Mathematica. 



III. SHORT AND LONG-TIME DIFFUSION 
COEFFICIENT 



In the following we evaluate the three contributions to the 
average velocity of the sphere, i.e., f3(D ss ) ■ F ext , (v I s ) F and 
~ fromEqs. {14} with (-} F = (-) F (t -> oo). Since 



3r\F 



(D ss ) is explicitly multiplied by Y ext , in the linear response 
regime the coefficient D ss in the first contribution has to be 
averaged with respect to the equilibrium distribution P eq for 
all times t. Therefore the short time diffusion coefficient D s s 
of the sphere equals the first contribution to the long time dif- 
fusion coefficient 

oo 

g = 1 + ^LpJdrr 2 [{D: s {t) 1) + 2(D?,(r) - 1)] 



= 1 - a s p, 



(21) 



with 



4 



3a, = 5R 3 5R g R 3 2Rj + 5R 2 g R 3 3 -^f - 12R g B* - 10R 3 g R 3 
4tt 2 2R R 2 R 3 

6Q9RI + 7500R g Rl - 4£>mR 2 g R\ - 875R 4 g R 3 h -2079R g Rl - ll250R 2 g R 6 h + 56QQR 3 g Rl + 175R*Rl 

140i? 4 140i? 5 
-36RI + 882R 2 g R 7 h + 3Q00g|fl| - 1225^ _ 36R g Rj - 294R 3 g R 7 h - 750R i g R e h + 245g|g| 

56i? 6 56i? 7 ' ( ^ 



For evaluating the interaction velocity defined in Eq. ( fT5l ) in 
the stationary limit, we note that the gradient of the potential 
V is only nonzero at r = R and points in direction r. This 
leads to 



Air 

«) F = f3D° sY R 2 A(R)(Dl s (R) 
= -cuf3D°F ext . 



D 



(23) 



With a partial integration with respect to r the stationary 
Brownian velocity defined in Eq. ( fToT ) is given by 



<vfT 



f3D° s ^-p J drA(r) 







-r'(Dl s (r)-D r sp (r)) 



-2r(D s ss (r)-D%(r)) 
= -a B rpD° s F ext . 



(24) 



Combining the three contributions, the long time diffusion co- 
efficient of the sphere is finally given by 



Di 



D° [1 - (a s + ai + a Br )p] = D° s (1 - a lP ) . (25) 



Eqs. ( TZTT i and ( fZ5b can also be stated as functions of mass 
density c (which is often used in the experimental literature) 
of polymer coils rather than the number density p by using 
c = p-fif-, with Avogadro's number Na and the molecular 
mass of the polymer coils M. The coefficients a' s and a[ re- 
lating the mass density to the short and long time diffusion 



constants are then given by a' s 
spectively. 



■ ^W and a 'i 



at 



Na 
M 



IV. RESULTS 



The coefficients a s and a; for the short and long time coef- 
ficient are functions of R s , R g and Rh and they can be writ- 
ten as the product of R 3 and a function that depends only 
on ( = R g /R s and £ = R g /Rh. For hard spheres £ = 1. 
For large polymers in good solvent conditions, £ approaches a 
universal value of £ ps 1.6, see ll36tl and references therein. 
In order to be able to compare our results to the case of 
hard spheres J5|, we introduce the polymer packing fraction 



^-R 3 g p and define 



a s /t(R s ,R g , R h )p = 5 a /i(C, £) 4> P - 
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FIG. 3: ai and a s (solid blue and dashed red lines, respectively) 
defined via D l J s / D a s = 1 — <ii/ s <j> v as functions of C, = R g /R s 
for different values of the ratio £ = R g /Rh- The value of £ for 
neighboring curves for the lowest six curves differs by 0.1. Squares 
and circles indicate the results for the long and short time diffusion 
constant of hard spheres (£ = 1), respectively, according to Ref. 0]. 



Fig. [3]shows a s (£, £) and a;(C, £) as function of C for different 
values of £ between one and two. 

For £ — > 0, i.e., for tracer particles large as compared to 
the polymer coils (this limit is often referred to as the col- 
loid limit), the polymer solution as seen from the colloid be- 
haves like a continuum, the probability distribution P°°(r) 
approaches P eq (r) and both, a s and a; converge to the con- 
tinuum result 



lim a s 



= lim ai = --r. 

C^o 2£ 3 



(27) 



This leads to the generalized (sometimes also called effective) 
Stokes-Einstein relation for the diffusion coefficients 

k B T 



lim D\ = lim Dt 



(28) 



67ri? s ?7 00 

with the Einstein result for the zero-shear high-frequency lim 
iting viscosity of the polymer solution 



l-i^-Rl 



2 3 



hP 



= m 



l 



5 47T 

2 ~3~ 



R 3 hP +0( P '). (29) 



?7o is the Newtonian viscosity of the (polymer free) solvent. 
Note that the next term in a s /i is of 0(C), while the difference 



(26) between a s and a; is of 0(C ). 
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FIG. 4: ai and a s (solid blue and dashed red line, respectively) for 
large polymers in good solvent conditions (£ = R g /Rh = 1.6) as 
function of f, Also shown is the result for the generalized Stokes- 
Einstein relation Eq. J28b which is independent of £ = R g /R s . In 
this approximation short and long time diffusion constant are equal. 



For C — * oo, i.e., in the so called protein limit in which 
the tracer particle is small as compared to the polymer coils, 
a s — > and the short time diffusion coefficient approaches 
the value in the pure solvent D a s from below. The long time 
diffusion coefficient reaches a finite value which is smaller 
than D°- 



lim ai = \. 

C— too Z 



(30) 



It is important to note that our model is limited to the case 
were the colloid does not enter the polymer, i.e., the colloid 
must remain larger then the "mesh-size" of the polymer. In 
the protein limit, the diffusion coefficients are independent of 
hydrodynamic interactions since the small tracer does not per- 
turb the solvent significantly. Both short and long time coeffi- 
cient approach the value of the corresponding calculation ne- 
glecting hydrodynamic interactions. This simplifies Eq. (jTTJ 
significantly and we get the analytic results 54 noHI) = and 

-(aoHI) = (1 + 1/C) 3 /(2 + 2e/C)iWith£l (noHI) ^ 1 ^ ( ^ ^ 

For a very long polymer in good solvent conditions £ w 1.6. 
The corresponding values of a; and a s as functions of ( are 
shown in Fig. [4] In contrast to the short time coefficient a s 
which decreases monotonically as a function of £ the long 
time coefficient 5/(£,£ = 1-6) has a maximum at ( w (. 
That means that for Rh < R s (£ < £), larger spheres 
are less hindered in their motion by the polymer coils than 
smaller spheres, while the situation is reverse for Rh > R s 
(C ^ O- This is in agreement with experiments as demon- 
strated in Sec. IV A| The variation of 5/ over all £ is neverthe- 
less rather weak (about ±15%), which means that the general- 
ized Stokes-Einstein relation is an acceptable approximation 
for all values of £. However, this is only the case for £ rs 1.6, 
for which the limits of hi for £ — > and for £ — > oo are not too 
different. In contrast to the long time diffusion coefficient, the 
short time diffusion coefficient varies strongly as a function 



of £, such that the generalized Stokes-Einstein relation for fi s 
holds only in the limit £ — > 0. The maximum of a; (£, £ = 1.6) 
at £ = £ is the result of two competing effects. With decreas- 
ing R s (increasing £) the solvent flow field generated by the 
moving tracer is weaker and as a consequence the short time 
diffusivity increases, i.e., a s and, according to Eq. dZST l, 5/ 
decrease. On the other hand, the distribution of bath parti- 
cles around the tracer gets more disturbed since the weaker 
flow field cannot transport the bath particles around the tracer 
such that these accumulate infront of the tracer, which reduces 
the tracer mobility ITil [l5ll . For very large £ the decreasing 
short time coefficient 5 S dominates but at intermediate £ the 
accumulation of bath particles leads to a local maximum of 
ai. This mechanism only leads to a local maximum of a; for 
intermediate values of £. 

The long time diffusion coefficient for a tracer particle in a 
suspension of equal hard spheres is expected to be 5/(1, 1) = 
2.10, see, e.g., 112111 . Our theory yields a slightly larger value 
5j(l, 1) = 2.14 since we neglect lubrication forces at small 
particle distances (which are less important if the interaction 
radius is larger than the hydrodynamic radius, i.e., for £ > 1). 
For this reason a; for £ = 1, while still beeing monotonic, 
shows the onset of a local maximum in our theory (see Figj3l), 
which is in contrast to the results obtained in Refs. ifsL l37ll . 
For polymers with the same interaction radius as the tracer 
sphere we find Sj(l, 1.6) = 0.66. Therefore a tracer particle 
is much less hindered in its motion by a suspension of equal 
sized polymers than by a suspension of equal sized spheres. 
Because the polymer hydrodynamic radius is smaller than the 
hydrodynamic radius of the hard spheres, the tracer and the 
polymer interact less strongly via the solvent. 

In order to test the accuracy of our results obtained with 
hydrodynamic tensors up to order 0(r~ n ), we repeated the 
calculation with hydrodynamic tensors of the next lower or- 
der (0(r~ 9 )). The relative deviation A of the two results are 
A < 8% (£ = 1), A < 4% (£ = 1.25), A < 2% (£ = 1.6) 
and A < 1.1% (£ = 2) for both long and short time results. A 
vanishes for both ( — > and £ — > oo. As expected, the order 
of the expansion is less critical for £ > 1, i.e., for small Rh, 
and taking into account even higher orders should not change 
the results for di/ s significantly. Calculations with lower or- 
der approximations to the diffusion tensors (0(r~ 7 ) or less) 
do not yield the correct generalized Stokes-Einstein relation, 
Eq. <|28]). 



V. COMPARISON TO EXPERIMENTS 

A. Stretched exponential and scaling exponents 

Experimental values of the long time diffusion coefficient 
Ir4lll of tracer spheres in polymer solutions are often described 
empirically by a stretched exponential 



DO 



-Cc v M~<R° 



(31) 



with a dimensional constant C. It has been noticed that the 
form OH has unphysical limits for both 5 < and S > 



6 



for R s — > oo and R s — > 0, because one expects a finite value 
Di/D^ 7^ 1 for any i? s . Rescaled versions have been sug- 
gested I25I.I29I1 . where only the difference between the limits 
R s — * 00 and i? s — > are described by a stretched expo- 
nential. This rescaling makes a direct comparison of 5 to ex- 
periments difficult, since experimental values are usually ex- 
tracted from un-rescaled data. Apart from this, the general 
experimental findings for the protein limit (£ ^> 1) are S > 
(see Hzl Table 1]), i.e., D l s /D° decreases with R s . In the 
colloid limit (C < 1), D l JD° was found to be almost inde- 
pendent of R s with a negative S f=a —0.1 U. The maximum 
of di at £ rs 1 (see Fig. |4]i is therefore in qualitative agree- 
ment with experimental findings. Since a; is non-monotonic 
as a function of £ in our calculation, we conclude that D l s /D° s 
cannot be described by a stretched exponential in R s over 
the whole range of size ratios. This is in disagreement with 
Ref. i29ll . in which a universal value 5 — 0.77 and a mono- 
tonic behavior was found. This difference might be due to the 
fact that our prediction is valid in the dilute limit, while the 
model in l29ll is based on a depletion layer which is more pro- 
nounced at higher densities. We also emphasize that the short 
time diffusion coefficient decreases monotonically with R s . 
In Ref. l29ll . the distortion of P which gives rise to the differ- 
ence between short and long time diffusion was not taken into 
account. 

Let us turn to the other exponents in Eq. ( TJTl ). In the dilute 
limit, i.e., for small c — p M/Na, the exponent in Eq. OTb is 
small and D l s /D° s w 1. In this limit, for which our model is 
made, we get 

r>i 

j± = 1 - C c v Mi R? s + O {(C c v M~< R s s ) 2 } . (32) 

In terms of mass density c, our result for the long term diffu- 
sion constant in Eq. ( 125b reads 

4 = 1- ^~a lR lc^. (33) 

D Q 3 1 9 M 

As illustrated in Fig. [4] the dependence of 5; on £, i.e., on the 
ratio of the polymer radius of gyration to the size of the tracer 
particle is rather weak. We therefore neglect this dependence. 
Using the scaling of the polymer size with its mass, R g oc 
M"» JH], we find by comparing Eqs. d32b and d33l 

7 = 3v g - 1. (34) 

In a good solvent, self avoiding walk statistics for a Gaussian 
chain lead to v g = 0.588 for the average size of the polymer 
ll38tl . Using this value, we find 

7 = 0.76, (35) 

with good agreement to the experimental value 7 = 0.8 found 
in Ref. |9|]. Note that the generalized Stokes-Einstein rela- 
tion also leads to the expression in Eq. d34l ) for 7. The ex- 
ponent v of the concentration c is in our linear theory equal 
to unity by construction. In summary, all exponents compare 
well to the experimental findings, see Tab. Q] For polymers at 
0-conditions with v g = 0.5 we predict 7 = 0.5. 



V 


7 


<5 


Rg Rs 






Experiment 0.6 . . . 1.0 


0.8 ± 0.1 


-0.1... 0.0 


Eq. {33} 1 


0.76 


weak dependence 


Rg ^ Rs 






Experiment [27] 0.5 . . . 1.0 




0.69... 1.0 


Eq. (33J 1 


0.76 


S > 



TABLE I: The exponents from Eq. t3U as measured in experiments 
compared to our theoretical prediction given by Eq. d33t . 



The phenomenological law OTb is best valid at semi-dilute 
polymer concentrations as stated in Ref. @|. It is in fact non- 
analytic in c for v ^ 1. Despite this, the experimental data 
summarized in ^ also covers the dilute regime and shows 
similar behavior there, see Sec. IV Bl However, we found only 
very few experiments fioll focusing on the dilute regime. 



B. Quantitative comparison 

Our model does neither include long ranged forces between 
the tracer particle and the polymer coils nor adsorption of the 
polymer on the tracer. This situation is realized in Ref. |8y: 
poly(ethylene oxide) (PEO) is a neutral polymer, so electro- 
static interactions between the sphere and the polymers are 
absent, and polymer adsorption on the sphere is suppressed 
by a surfactant. 

Neither the radius of gyration R g nor the hydrodynamic ra- 
dius Rh of the PEO polymers were measured in the experi- 
ments and we calculate them according to Ref. l39ll via 

R g = 0.0215A/ a583 nm, (36) 
R h = 0.0145M°- 5n nm. (37) 

Fig. [5] compares the diffusion constants measured by light 
scattering Q Fig. 3c] with our results, which are by construc- 
tion linear in the polymer mass density c. The tracer sphere 
(polystyrene) has a radius of R s = 322 nm and the poly- 
mers have molecular masses of 18500, 10 5 and 3 ■ 10 5 amu. 
The experimental data for the smallest polymer with M = 
7500 amu in iH Fig. 3c] were not useful at small concen- 
trations due to large scatter. According to |@] the diffu- 
sion constant of the tracer particle in the pure solvent was 
= 7.44 • 10~ 9 cm 2 s _1 . Given these values, there is no fit 
parameter for the initial slopes in Fig. [5] The overlap concen- 
tration c* = 3M/(4nR 3 g N A ) is 25.5, 7.1 and 3.2 g/L for the 
three polymers, respectively. Our theory, which is linear in the 
polymer concentration, is valid only for c <C c*. In Fig. [5] the 
experimental points start to deviate from the straight line for 
smaller and smaller c as M increases. For M = 3 • 10 5 amu, 
they deviate considerably for c > 1 g/L. The agreement for 
concentrations much smaller than c* is good in all three cases. 
Note that according to our theory the short time diffusion con- 
stant Dg is almost identical to D l s since the polymer coils are 
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FIG. 5: Diffusion coefficient D s of a polystyrene sphere of radius 
R s = 322 nm in a solution of PEO polymers with molecular masses 
of 18500 amu (squares), 10 amu (circles) and 3 • 10 5 amu (trian- 
gles) from Ref. [g]. The solid, dashed, and dash-dotted lines, re- 
spectively, indicate the theoretical predictions for D l s according to 
Eq. j25\ . There is no adjustable parameter. 

small as compared to the tracer particle, i.e., ( <C 1, for the 
cases shown. 

Fig. [6] shows the normalized diffusion coefficient for two 
different sphere sizes (322 and 51.7 nm with D a s = 4.64 ■ 
10~ 8 cm 2 s _1 JH]) in a solution with a polymer mass of 
3 ■ 10 5 amu. For small polymer concentrations the normal- 
ized diffusion coefficient depends only weakly on the sphere 
size. For the smaller sphere, the ratio of R g to R s is given 
by C = 0.65, such that a continuum theory is not expected to 
hold. For this value, long and short time coefficient should 
differ appreciably but the experimental values lie between our 
predictions for the short and long time diffusion constants. 
From Ref. |8[] it is not clear weather the experiments probe 
the long or the short time coefficient. Note that at higher con- 
centrations, the larger sphere is less hindered in the diffusion 
by the polymers, which is consistent with a negative value of 
S. 



affine way on the polymer density. While the short time diffu- 
sion coefficient decreases monotonically as a function of the 
size of the tracer particle, polymers with hydrodynamic radius 




0.5 1 1.5 2 2.5 3 3.5 4 
c [g/L] 

FIG. 6: Normalized diffusion coefficients D s of spheres with radius 
R s = 322 nm and 51.7 nm (circle and square, respectively) in a 
solution of PEO with a molecular mass of 3 ■ 10 s amu from Ref. @]. 
Dashed and dash-dotted lines are theoretical predictions for D l 3 and 
for R 3 — 51.7 nm in first order in polymer concentration from 
Eq. {25) and Eq. {2D, respectively. For R s = 322 nm, D l s and D% 
are almost identical (corresponding to £ = 0.1, see Fig. |4j and we 
only plot D l s (solid line). 

comparable to the tracer size seem to be most efficient in de- 
creasing the long time diffusion coefficient. 

Our current model is limited to low polymer densities, 
but direct polymer-polymer interactions can be taken into ac- 
count in the framework of a dynamic density functional the- 
ory lfl2tl . Hydrodynamic interactions between the polymers 
and the tracer sphere can be taken into account in the same 
way as in this paper |fl4tl and recently dynamic density func- 
tional theory has been extended in order to take into account 
hydrodynamic interactions among the polymer coils l40ll . 
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VI. SUMMARY 

We developed expressions for the short time and the long 
time diffusion coefficient of a tracer sphere in a dilute solution 
of particles with a hydrodynamic radius which can be smaller 
than the hard sphere radius for the interaction with the tracer 
particle, e.g., polymers. Solvent mediated hydrodynamic in- 
teractions are taken into account up to 1 1th order in reciprocal 
distance, therefore neglecting lubrication forces at close dis- 
tances. Calculating the diffusion coefficients is reduced to the 
solution of a system of 1 1 algebraic equations. 

The results are in good agreement with experiments in the 
dilute regime, in which the diffusion constant depends in an 
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APPENDIX A 

The coefficients D^(r), a £ {r , 8} and i,j £ {s,p}, of 
the diffusivity matrix for a polymer (hydrodynamic radius Rh) 
and a spherical tracer particle of hydrodynamic radius R s read 
up to order r" 11 (Ref. JH) 
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D r ss (r) = 1 - %S%R. (J)' + - 2RlR?) (i)' - (3^ - 7M£f& + 5i«) (;)' 

- ji^-R s (9.R£ - l20R 2 s R 2 h + 375R 3 s R h + 70R 4 S ) Q + O ^ Q , (Al) 
D e ss (r) = 1 - ^R 5 h R a (I) - Ir s (9Rl - 9R 2 s R 5 h + 10i«) (J) 

-^i? s (6i?£ - 18R 2 s R 7 h + 35i?X) Q + ((J) V (A2) 
^p(r) = 1 - ^R h Rl (J)' + (^f^ - JUfeitf) (1)' - ^ (3i?I - 22i« + 5i«) (i)' 

~R h Rl (70R 4 h + 375R s R 3 h - 120R 2 R 2 h + 9i£) (^j + O ^ Q , (A3) 
D e pp (r) = 1 - ^R h Rl (jj - ^R h (9R 7 S - 9R 2 h R b s + WRfM) (J) 

~R h (6R 9 S - l8RlR 7 s + 35R 4 h R 5 s ) (J^j +0^ Q ^ , (A4) 

W = ^-\Rs (Rl + Ri) (J)' + ™RlK Q) 7 - f i« + Rl) (J) 9 

+ ^ (104 - lh\R 2 s Rl + 10i?X) (J) + O ^ ^) j , (A5) 

= ^ + J^ + ^)(J) 3 + 4^8o^^ (A6) 
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